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We present two universal models of quantum computation with a time-independent, frustration- 
free Hamiltonian. The first construction uses 3-local (qubit) projectors, and the second one requires 
only 2-local qubit-qutrit projectors. We build on Feynman's Hamiltonian computer idea [3 and 
use a railroad- switch type clock register. The resources required to simulate a quantum circuit with 
L gates in this model are 0{L) small-dimensional quantum systems (qubits or qutrits), a time- 
. . . independent Hamiltonian composed of 0(L) local, constant norm, projector terms, the possibility 

0^ \ to prepare computational basis product states, a running time 0{L\og^ L), and the possibility to 

. measure a few qubits in the computational basis. Our models also give a simplified proof of the 

' universality of 3-local Adiabatic Quantum Computation. 
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I. INTRODUCTION 



Since the introduction of the standard circuit model [H, many new models of quantum computation have 
r~i been developed. Inspired by Feynman's ideas about how to use quantum systems to simulate reversible computation 
'"5 \ 01 9- class of 'analog' Hamiltonian quantum computing models computing is based on time evolution with time- 
independent or slowly changing Hamiltonians. The basic idea is to encode the progress of a quantum computation U 
^ (a quantum circuit composed of a sequence of local unitaries) with L gates 

§■ U = Ul---U2Ui (1) 



on an initial state ^ into a history state of a system with two registers TLwork ® 'Hdock as 



= -=^^(C/i...f/o|^o)„)®|t),. (2) 

(N ; 

, First, one can now look for a simple Hamiltonian with ^ as its ground state. Finding the ground state of such a 
00 Hamiltonian, measuring its clock register and obtaining the result L would then allow us to obtain \(Pl), the result 
of the quantum circuit U applied to the initial state \ipo) of the work register. This approach produced many results 
about the quantum complexity (in particular, QMA-completeness) of finding the ground state properties of simple 
S< , Hamiltonians [1, i, 0, [11 E 111 [ll ■ 

J> , A second direction of investigation in these models [l^, [l^ [3 is to look at the Hamiltonian dynamics instead 
of the static ground state properties. The Hamiltonians in these Hamiltonian Quantum Computing (HQC) models 
generate Schrodinger time evolution which brings an easily-prepared initial state towards the history state ^ or its 
• variants. As we show in Section |TT1 one can extract the result of the computation U \(Pq) from the final state with 
high probability. This results in BQP universality - any problem which is efficiently solvable on a quantum computer 
can be solved in the HQC model as well. Another example utilizing the preparation of the state ^ is the result of 
Aharonov et. al. [l9l|, showing that Adiabatic Quantum Computation ^ is polynomially equivalent to the standard 
circuit model. 

In this paper, we take the second point of view, focusing on the universality of the dynamical properties of a class of 
Hamiltonians. The desirable properties for such a HQC model are low locality (a Hamiltonian composed of terms which 
act nontrivially on a few qubits at a time), a simple geometry of interactions, low dimensionality of particles involved 
(qubits or qudits), time invariance (or slow change with time), translational-invariance and a large eigenvalue gap 
above the ground state energy (or above the computational subspace) to protect the computation against decoherence. 
We achieve this last property by using a frustration-free Hamiltonian - a sum of positive semidefinite terms, whose 
ground state is the common ground state of all the terms. 
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FIG. 1: a) Computation on a "line" of states \tiit)- b) The line corresponding to a circuit U padded with extra identity gates. 
The region where the computation is done is marked, c) Computation on a cycle of states |i/)t) coming from two copies of the 
computation on a line, d) The cycle corresponding to a padded circuit U . 



Our goal is to find a universal HQC with a low locality of interactions, using a Hamiltonian that is time-independent 
and frustration- free. Feynman's original construction requires 4- local terms (see Section |lT] for details). We equip our 
quantum computer with a different, railroad- switch clock register. This results in an effective translation from the 
circuit model to a Hamiltonian (continuous-time) model of quantum computation with a simple, frustration-free, 
time- independent, now only 3-local Hamiltonian. Note, that to show the universality of our model, we do not need 
to use the adiabatic theorem, because our Hamiltonian is frustration- free (see Section llHBp . 

Recently, Mizel et. al. [l^ took an alternative approach with a construction using physical particles jumping down 
a grid, utilizing only 2-particle (but 4-site) interactions. For comparison, when wc look at their construction encoded 
using qubits, the interactions become 4-local. 

After finishing the present article, we found out that Q also contains a 3-local version of this computer, using a 
so-called switch primitive, unknown to the wider community. The details of the dynamics of Feynman's model with 
such a clock were then analyzed by DeFalco et al. in [l^, [2l| , later focusing on the entropy of the clock register in 
[13, H^. Nevertheless, there are two things that make our re-discovery different from previous work. First, we in 
the 3-local construction we focus on choosing a random time for the final measurement, which requires less detailed 
analysis than the one in [20| . We also look at the problem from a modern quantum information perspective, and 
make connections of our work to other universality models such as adiabatic quantum computation. Second, we take 
a step further and make the clock construction a two- local one (for a qubit and a qutrit). This paves the way for a 
possible future new proof of universality of 2-local Hamiltonians without the requirement of large norm penalty terms 
dictated by perturbation theory gadgets as in |10|] . 

The paper is organized as follows. First, in Section |lT] we review the continuous-time quantum computing model 
based on Feynman's idea and the usual ways of encoding the clock register. We present our main innovation, the clock 
register with railroad-switch gadgets in Section [Till and bound the required running time of the resulting 3-local HQC 
model in Section llH Al Next, we discuss how our result relates to the universality of Adiabatic Quantum Computation 
in Section IIIIBI and analyze the gap protecting the computational subspace in Section IffI CI Section IIVI contains our 
second result, a universal HQC model involving only 2-local qubit-qutrit interactions. Finally, in Appendix |^ we 
prove a Lemma about the mixing of a continuous-time quantum walk on a cycle, utilized in Section IHI Al for the 
analysis of the required running time of our model. 



II. FEYNMAN'S COMPUTER AND CLOCK CONSTRUCTIONS 



Already in 1985, Feynman [7[ investigated the possibility of simulating a sequence of reversible classical gates in 
a quantum-mechanical system. The construction he found applies equally well to implementing a quantum circuit, 
which is also reversible. Let U = Ul ■ ■ - Ui ^ be a sequence of L unitary transformations Ut on n qubits. Consider 
a quantum system with two registers, Ti-work Ti-dock: where the first of these holds n work qubits qi, . . . ,qn- The 



3 



role of the clock register is to hold a pointer - indicating how far the computation has progressed. A simple time- 
independent Hamiltonian then facilitates the evaluation of the circuit U on the work register of this system. Consider 
the Hamiltonian 



(3) 



where Ut acts on the work register. The first type of transition given by this Hamiltonian increases the clock register 
from \t — 1)^ to \t)^, while applying the unitary Ut to the corresponding work qubits. Similarly, the other transition 
decreases the clock register and uncomputes the unitary Ut- Fixing an initial state \fo)w of the work qubits and 
taking t = 0, . . . , L, we can define a basis 



|V't) = (C/t...C/2C/i|^o)^)®|t), 



(4) 



of a particular subspace Ti-ip^ of the whole Hilbert space Ti-work ® 'Hdock- Within each such subspace, Feynman's 
Hamiltonian has a very simple form: 



Hp 



1 

1 1 
1 



1 

1 



(5) 



It is a Hamiltonian of a quantum walk on the "line" of states \ipt), as each state {ipt) is connected by a transition in Hp 
only to its direct neighbors (see Figure [1]). Time evolution with Hp according to the Schrodinger equation thus brings 
the initial state |V'o) = Ivo)™ \0)c iiito some superposition of the states IV't) Q- After some time r, let us measure 
the clock register. When we obtain L, the state of the system after the measurement becomes = I'Pl)^ ® \-^)c- 
The work register then contains the desired output of the circuit U . One might argue that this doesn't happen too 
often. However, we can boost the probability of obtaining by slightly modifying our system. First, we pad the 
circuit with 2L extra identity operations 



U' = 



■3LJ13L-1 



(6) 



and, accordingly, expand the clock register to hold states |0)^ , . . . , |3L)^ (see Figure [TJd). All the states Ift)^ of the 
work register with t > L are now the same - equal to the desired output state \'Pl)w We now let the system evolve 
from lipo) for a time r chosen uniformly at random between and 0{Llog L). Because the quantum walk on a line 
is rapidly mixing, the probability to measure a state with the clock register \t > L)^ at time t is close to |, as we 
show in Appendix lAl This yields the desired state \(Pl)w of the work qubits. If a measurement of the clock register 
gives us < < L, we repeat the experiment. 

Note also that we can turn the Feynman Hamiltonian ([3]) into a sum of projector terms. To do this, we change the 
sign of Hp and add terms of the type \t — 1) (t — 1|^ and \t) {t\^ as 



(7) 



Just as in the Hamiltonian H^°-' restricted to a subspace spanned by the basis ^ has a simple matrix form 



1 -1 

-1 2 -1 

-1 2 -1 

-1 ■■. 



2 -1 
-1 1 



(8) 
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Observe that as the rows sum to zero, for any initial state of the work qubits \<fo)w^ corresponding history state 
-7== li/jj) ^ is the ground state of H^°-' , with energy equal to zero. Again, H^°-' is a Hamiltonian of a 

quantum walk on a line, although now with a small modification at the endpoints and a shifted spectrum: 

Hr'^l-lHF-lm{0\^ + \L){Ll}. (9) 

This is the form of the Hamiltonian which Kitaev used in (along with additional terms) to prove QMA-completeness 
of the Local Hamiltonian problem. Finally, note that we can get rid of the cndpoint terms by wrapping the computation 
around on a circle as in Figure [TJ;. 



A. Implementing the Clock Register 

So far, we haven't said anything about how the clock register is actually implemented. Its states should indicate 
where in the computation we are, encoding the L + I different states \t)^. The first naive approach is to use a binary 
encoding of t which would require log(L + 1) qubits. However, the transition operators \t) (t — 1|^ in the Feynman 
Hamiltonian would then be highly nonlocal. It is much better to use the idea Feynman originally had in mind. Imagine 
the clock register as the state space of a pointer particle hopping on a line 0, . . . , L, with the state \t)^ corresponding 
to the pointer positioned at t. This pulse clock uses L + 1 clock qubits cq, . . . , cl, representing the states \t)^ as 



|0), = 


|1000. 


..00), 


(10) 


|l>c = 


|0100. 


..00), 




|2). = 


|0010. 


..00), 






|0000. 


..01), 





with a single clock qubit (indicating the pointer) in the state 1 1). What is the reason to take such a generous approach, 
using only L + 1 of the 2^+^ states available in the state space of the pulse clock register? For the pulse clock, the 
transition operators in ([3]) acting on the clock register can be implemented just 2-locally, acting nontrivially on only 
two clock qubits and as an identity on the rest of the system: 

K)(i-l|^ = I® |01)(10|^^_^^^^®I. (11) 

Two-qubit unitaries Ut are necessary for universal quantum computation, so the Hamiltonian ([3]) with the pulse clock 
encoding contains 4-local terms. Our main result, given in Section lllli uses a pulse-type clock and implements dU 
just 3-locally. Note though, that there are two problems with the pulse clock. First, it is not entirely simple to locally 
check the clock register for correctness (ruling out the states with two neighboring I's, but also with two far-away 
I's). Second, the pulse clock requires initialization - ensuring that the clock register has a pointer (a spin in the state 
|1)) at all. This can't be done by using only local positive semidefinite terms in the Hamiltonian, which makes the 
pulse clock unsuitable for obtaining QMA^ -completeness results such as [l2l. [l3. [25| . 

Another way of encoding the states \t)^ is a domain wall clock, using a unary representation of |t)^, introduced by 
Kitaev [Sj for his Local Hamiltonian problem. The progression of states (on L + 2 clock qubits) is 



|0), = 


|1000., 


. .000) , 


(12 


|l>c = 


|1100., 


. .000) , 




|2)e = 


|1110., 


. .000) , 






Illll . , 


..110), 





where the "clock time" (or alternatively, pointer position) t is determined by the position of the domain wall 10. A 
useful property of the domain wall clock is that we can check whether the state of the clock register is in one of the 
states \t)^ (|12l) by 2-local operators. What suffices is to ensure that the sequence 01 never occurs. Moreover, we 
can ensure that there is a single domain wall 10 in the clock register by imposing that the first qubit is fixed to 1 
and that the last qubit is 0. This solves the initialization problem of the pulse clock (where we can't rule out the 
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FIG. 2: The pulse clock is a line (a single track) of L + 1 clock qubits with a single active qubit (the train) in the state |l). 
When the clock qubit ct is active, the pulse clock encodes the state \t)^. 

state |00 • • • 00)^ by simultaneously satisfiable local positive-semidefinite operators). With the domain wall clock, the 
transition operators in Hp ^ arc 3-local, as we need 

\t) {t-ll = l® |110) (100|^^_^^^^_^^^^ ® I. (13) 

Recall that we want to increase the clock register while appling two-qubit unitaries. Feynman's Hamiltonian with 
the domain wall clock is thus 5-local. If one uses transition operators with smaller locality (as in [H, [l3l), transitions 
to illegal clock states (states of the clock register other than the prescribed \t) ^ will occur in the Hamiltonian time 
evolution with Hp- This needs to be dealt with by imposing large energy penalty terms in the Hamiltonian. Moreover, 
the Hamiltonian will no longer be frustration-free. For this reason, the domain wall clock does not fit our present 
purpose: finding a simple, low-locality, frustration-free Hamiltonian with a small norm, allowing us to do universal 
quantum computation by the time-evolution it generates. 

Another approach is to combine the two clock constructions above by interspersing a domain wall clock lattice with 
a pulse clock lattice. In this fashion, Nagaj and Mozes [l2| and Eldar et al. [13] showed QMA-completeness of the 
3-local Hamiltonian problem [2BI with a much larger promise gap. 

Going beyond the linear structure of the pulse/domain wall clocks, Eldar et al. (T^ have recently come up with 
a novel idea. They use a progression of states of the clock register that can go from \t — 1)^ to \t)^ by two different 
routes, depending on the state of one of the work qubits. This results in their proof that the Quantum 2-SAT problem 
[25! for a cinquit-qutrit pair (each Q-2-SAT condition acts on one 5-dimensional and one 3-dimensional quantum 
system) is QMA^-complete [l3|. In the next section, wc take the nonlinear route as well, construct a novel, pulse-type 
clock and make the Feynman Hamiltonian 3-local (as opposed to 4- local with the usual pulse clock) . What makes our 
construction different from is a bit-flip symmetry which allows us to easily analyze the system's dynamics and 
prove its universality. 

III. THE RAILROAD SWITCH: A 3-LOCAL GADGET 

For a given quantum circuit J7, our goal is to construct a Hamiltonian which facilitates time evolution of a particular 
initial state towards a state very much like ([2]), encoding the result of U . Specifically, we focus on the clock register 
encoding of the Feynman Hamiltonian ([3]) . We base our construction on a pulse clock pO)) . We can view the progression 
of the pointer 1 in the clock states |t)^ as a train moving along a single track as in Figure^] As it goes between stations 
Ct_i to Ct, one-qubit or two-qubit unitary operations Ut are applied to the work qubits. For any single qubit gate, the 
operator Ut ® \t) {t — 1|^ is already 3-local pTjl . Besides single-qubit gates, all that is necessary for universality is the 
ability to apply the 2-qubit CNOT gate. As we have seen in the previous Section, 4-local interactions are required 
for this in the standard Feynman Hamiltonian construction. To do it with only 3-local interactions, we have to add a 
twist to the pulse clock, going beyond its linear structure. For each CNOT gate, instead of viewing the clock register 
as a train on a single track as in Figure [H we introduce a railroad switch gadget with two tracks, depicted in Figure 
[3l The train will use the upper or lower track depending on the state of a 'train master' qi - one of the work qubits. 
Only on the upper track, the target work qubit q2 is flipped. The gadget thus facilitates the application of a CNOT 
gate on the work qubits qi and 92 ■ 

Our resulting Hamiltonian H^s is composed from positive semidefinite 3-local terms, so it is a Quantum 3-SAT 
Hamiltonian. We prove the universality of this 3-local Hamiltonian Computer model where we let an easily prepared 
starting state evolve with H^s for not too long a time, and also for an Adiabatic Quantum Computation (AQC) 
where we slowly change the Hamiltonian from a simple starting Hamiltonian Hb to H^s- We discuss this connection 
to AQC in detail in Section HITbI 

Let us present the details of the railroad switch. The original pulse clock has L + 1 qubits cq,. ..,cl. When a 
CNOT gate is due at step t (i.e. C/t =CNOT), we introduce four extra clock qubits (mi, 1*2, ^1,^2) between ct-i and ct 
as in Figure [21 Let us call the control work qubit qi, and the target qubit for the CNOT g2- If 9i is in the state |1), 
the active site (the train) in the clock register can only follow the upper track. For qi in the state |0), it must take 
the lower track. Furthermore, when the train moves from ui to U2, we flip the target qubit g2- The train then moves 
to the clock qubit Ct and the two tracks (computational paths) merge. This effectively applies a CNOT gate between 
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U2 





— c. 



X^0)(0|, 



FIG. 3: The railroad switch gadget for the apphcation of a CNOT gate between work qubits qi and 52, using only 3-local terms. 
Comparing to the pulse clock in Figure [2] there are four extra clock qubits inserted between ct_i and Ct, creating two tracks. 
The routing of the train (active site) from Ct-i (and ct) to the upper and lower tracks is conditioned on the state of the train 
master - a control qubit qi in the work register. Furthermore, when the train moves from ui to U2 on the upper track, the 
target work qubit 92 is flipped. 



qi and q2. The transition operators for the upper track are now all 3-local: 

® {\u,){ct^,\ + \ct^,){u,\), (14) 

< ® {\U2)(UI\ + \UI)(U2\), (15) 

|1>(1|,, ® {\ct){u2\ + \u2){ct\), (16) 

because the operators (ct|, etc. can be implemented for our pulse-type clock only 2-locally as 

|ui)(q_i| - |10)(01|^^_^„^, (17) 

\u,){u2\ = |10)(01|„^_„^, (18) 

\u2){ct\ = |10)(01L^^,^, (19) 



acting trivially on the rest of the clock qubits. The Hamiltonian for the upper track is then a sum of p4)) - (|T6l) . If we 
desire so, we have the option of making it a sum of projectors by adding a few terms as in The projector version 
of the Hamiltonian for the upper track is then a sum of three 3-local projectors: 

Hi'^ = \l){X^li\ct-i){ct^i\ + \u,){u^\~\u,){ct-i\-\ct^i){u,\) (20) 

+ ^ ® ( |wi> ("1 1 + 1^2) (U2 I ) - ^ ^ {\U2) {ui\ - \ui) {U2 I ) 

+ |1> ®\{\U2) {U2\ + \ct) {Ct\ - \U2) {Ct\ - \ct) {U2\). 

Note that when the clock register is in a proper clock state, i.e. one and only one clock qubit is in the state the 
projectors such as |ct) (ct| are easily implemented 1-locally as |1) (Ij^j- 
Analogously, the transition operators for the lower track are 

|0)(0|,^ ® (|/i)(c*_i|-f |ct_i)(/i|), (21) 

K ® {\l2){h\ + \h){l2\), (22) 

|0)(0|,^ ® {\ct){l2\ + \h){ct\). (23) 
Again, we can add a few terms and turn them into projectors, obtaining 

HfJ = \^)m,,®\{\h){h\ + \ct-i){ct-i\-\h){ct-i\-\ct-i){h\), (24) 

+ ^( l^l) ^ll + 1^2) (/2I- 1^2) (^ll- 1^1)^2!), 

+ \^){%,®\{\l2){l2\ + \ct){c,\-\ct){l2\-\l2){ct\). 

We now replace each term -ffp^o? ^o^' those t where Ut =CNOT with the railway switch gadget Hamiltonian 



rr(«) _ TT(t) I TT(t) /2c;^ 
"-switch — "up ^ "lo ' V^°7 
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obtaining 

^35 - E + E ^sltcH (26) 

t:l-qubit Ut CNOTs 



A. The Dynamics of our Model 

Let us now look at the time evolution generated by this Hamiltonian, focusing first on a single railway switch. 
Denote the work qubits on which the CNOT gate acts qi and q2 . Consider the action of on the state 



IVt-i) = \^t-l)^u ® |t - 1>, = (a |0>„^ \a)„,^_ +6 |1>„^ \(3)„^_ ) ® \ct-i), , (27) 




where \ ft-i) = a \^o...) +b is a general state of the work register and the clock register has the active site (train) 

at Ct-i- While the train is within the railroad switch gadget, the transitions given in pS)) dictate that the state of 
the system does not leave the subspace spanned by the following four states: 

IV't-i) = (a|(^o...) + ^'lvi...)) ® |ct-i),, (28) 
= a\^o...)®\h), + b\^i...)®\u,)^, (29) 

IV'E) - a|^o...)®|Z2)e + fc|^'i...>®l"2)e, (30) 
l^t) = {a\^o...) + b\^[„))®\ct),, (31) 

where \(p'i ) = aqf' \ipi...)- Observe that the work qubits of the state \^pt) have the CNOT gate applied to qi,q2, as 
we wanted, meaning that 

IVt) = CNOT,„,J^t_i)^ ® |t)^ . (32) 

This gives us a method of applying a 2-local CNOT gate while updating the clock register at the same time - via 
only 3-local interactions between the work and clock registers. 

Before introducing the railroad switch, the progression of states . . . , \tpt~i), |^t), |'0t+i), ■ • • was a "line", 

as each state in it was connected only to its two nearest neighbors by a transition in H . Now, with the transitions 
given by (pS)) . the modified progression of states . . . , |V'ti\), |'!/'t^\), |V't)j • ■ • again forms a "line". Let us call 

Tiiegai = span{\ipt)} (33) 

the subspace of Ti, spanned by this new progression of states, including those of the form (|29p and ([50)1 . Because we 
constructed it so, H^s (EH) does not induce transitions between the subspace TLiegai and Hj^^^;. The time evolution 
of an initial state in Hiegai is thus governed solely by the restriction of H^s to Hiegai- Moreover, in the new, updated 
basis |'0t) this restriction has the form ([5]) which we already know. It is a Hamiltonian for a quantum walk on a line of 
length L + 1 = poly{n), where n is the number of qubits the quantum circuit U acts on. The railroad switch requires 
three steps for each CNOT gate, so L = Li + 3Lcnot, where Li is the number of single qubit gates and Lcnot is 
the number of CNOTs in the circuit U. To avoid complications in the analysis of the required running time coming 
from the endpoints, we can change the 'line' of states \ipt) into a circle of length 2L = 2{Li + SLcnot) as in Figure 
[T];. First, we double the clock register as 

\coCi...Cl) » |coCl . . . Cl) (g) |C0'C1' . . . CL') , (34) 

where the Hamiltonian involving the second clock register has the same form (j26p . Second, we identify the endpoint 
qubits, i.e. cq = cqi and cl = cl'- This gives the clock register the geometry of a cycle, with a unique state 



io): 



icoo::>) (35) 



corresponding to t = 0. The active site (spin up) in the clock register can proceed towards cl both ways, as in Figure 
[it. The span of this new set of states \ipt) with the geometry of a circle then defines the subspace 'H'legai- When 
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we restrict the new Hamiltonian ((26)) including the terms involving the clock qubits Cf to the subspace 'H°^g^i and 
express it in the basis {IV"*)}! we get 



= i (I - i?o) , (36) 



where Bo is the adjacency matrix for a cycle of length 2L. The dynamics of this system is a quantum walk on a cycle, 
which we analyze in detail in Appendix]^ There we prove 

Lemma 2. Consider a continuous time quantum walk on a cycle of length L (divisible by 4), where the Hamiltonian 
is the negative of the adjacency matrix for the cycle. Let the system evolve for a time t < Tq chosen uniformly at 
random, starting in a position basis state |0). The probability to measure a position state \t) farther than L/6 from 

the starting point (the farther two thirds of the cycle) is then bounded from below as Po > ^ ~ jj^ — O |^ ^'°g _ 

Thus, when starting from the state concentrated on a single site of the line and letting the sytem evolve for a time 
chosen uniformly at random between zero and a number not larger than 0(Llog^ L), the probability to find the state 
farther than L/3 from the starting point is close to |. This corresponds to finding a state with a spin up on clock 
qubit Ct>L/3 (lower part of the circle) or Ct'>L/3 (upper part of the circle), depicted as the success region in Figure 
[TJi. In both cases, because the circuit U was padded with identity gates as in ([H]), the state of the work qubits is then 
Q and contains the output of the quantum circuit U . 

Therefore, the Hamiltonian (j26p allows us to simulate a quantum computation U with Li + Lcnot gates (here Li 
is the number of single qubit gates and Lcnot counts the CNOT's), as follows: 

1. Construct a pulse clock register with C = Li + Lcnot + 1 qubits. For each CNOT gate, add 4 extra qubits for 
the railroad switch gadget and modify the Hamiltonian as described by 



2. Pad the clock register with extra 2(Li + Lcnot) qubits and pad the circuit U with identity gates to make it 
three times longer. 

3. Double the clock register and join its endpoints. 

4. Initialize the system in the state 

m = \00...0}^^\0)l. (37) 



5. Let it evolve with H^g (given by (pS)) including terms for the whole clock cycle) for a time chosen uniformly at 
random between zero and r < O(Llog^L). 

6. Measure the clock register qubits in the success region (farther than L/3 from the initial site, see Figure [TJi). 
With probability close to |, we will find the active site (spin up) there. The work register now contains the 
output of the quantum circuit U. Restart otherwise. 

The Quantum 3-SAT problem ("Does a common zero-energy ground state of a sum of 3-local projectors exist?") 
is certainly NP-complete, but we still do not know whether it is also QM A ^ -complete [2^. Our railroad switch 
clock is not suitable for proofs of QMAi hardness of Quantum-fc-SAT, because it requires initialization (see Section 
HI)) . However, in this work we found that time- independent Quantum 3-SAT Hamiltonians can be universal in the 
following sense. We constructed H^g- a frustration-free, Quantum 3-SAT Hamiltonian (a sum of 0{L) 3-local projector 
terms) powerful enough to perform universal quantum computation in the Hamiltonian Quantum Computer model 
based on a quantum walk. When the number of qubits and gates involved in the circuit is L, what we need for our 
Quantum 3-SAT Hamiltonian computer are thus 0{L) qubits, a Hamiltonian with 0{L) projector terms with norm 
0(1), resulting in = 0{L), and a running time of order O(Llog^L). The rescaled required resources (time x 
energy) scale as r • \\H\\ — 0{L'^ log^ L). 

Furthermore, we show in Section Fill CI that the computation is protected by an energy gap scaling like 0{L~^) 
between the legal subspace Ti-^^g^i and the subspaces orthogonal to it, with the exception of the computationally dead 
subspace without a train (active site in the clock register). This energy barrier could be increased (at the cost of 
increasing the locality of interactions) by using error correcting codes like the ones developed for AQC by Jordan et 
al. [13. 

Although local (in the number of particles in each interaction term) , Feynman-like Hamiltonians are not geometri- 
cally local. Therefore, we would like to continue looking for Hamiltonians with geometrically local interactions (e.g. 
on a 2D grid), generating time evolutions towards a history state ^ or its variants, or usable in an AQC algorithm. 
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Oliveira and Terhal [ll| found such a 2-local Hamiltonian on a 2D grid. The downside of their construction is that 
the Hamiltonian doesn't keep the computational subspace invariant (it is also not is not frustration-free, with ground 
states which are history states). It thus doesn't generate evolution directly towards ([2|). Moreover, because it is 
based on perturbation gadgets, the eigenvalue gap of the final rescaled Hamiltonian itself is a rather large inverse 
polynomial in L. Its use for AQC is thus impractical. Two of the current ways of approaching this problem are then 
using higher-dimensional particles (qudits) or higher locality of interactions [18]. 



The model of computation we just described falls into a general category of Hamiltonian Quantum Computers 
(HQC) which initialize a system in a simple initial state, let it evolve for not too long a time, and end with a 
measurement of a few qubits (in the computational basis). Our HQC requires initializing the system in the state 
IV'o) (|37p . and then releasing it to evolve with the Hamiltonian H^s- We now look at the initialization in more detail 
and investigate how fast can we peform the 'release'. For the initialization, it is useful to think of ji/^o) as the unique 
ground state of the 2-local, initial Hamiltonian 



which is diagonal in the computational basis. While Hdockmit prefers the clock register in the state \t = 0)^, the term 
Hdatainit givcs an energy penalty to states whose work register is not properly initialized, i.e. all zero when the clock 
is in |co)j,. To let the HQC run, we turn off Hinu and turn on H^s- One way to do it is linearly as 



For our choice of H and a large enough Tmit, the conditions of the adiabatic theorem (see e.g. [28|) could hold. 
This model would then become an Adiabatic Quantum Computation (AQC) [5], where the system stays close to the 
instantaneous ground state of the time-dependent Hamiltonian H(t). At t = Tinu, the system would thus be close to 
the ground state of H^s - the history state ([2]) . 

Because of the energy gap, the AQC model has some protection against noise during the computation. The 
robustness and fault _tolerance of the AQC has been recently investigated by Childs et al. [1^ , Jordan et al. [131 , 
Lidar and Lloyd [30]. The challenge today is to find more new, entirely adiabatic algorithms, such as the adiabatic 
state preparation of Aharonov and Ta-Shma |31|]. 

In [19|, Aharonov et al. proved that Adiabatic Quantum Computation is polynomially equivalent to the standard 
quantum circuit model of quantum computation. This proof has two directions. First, as shown before by Farhi 
et. al. (23 | we can simulate time evolution with an AQC Hamiltonian on a standard quantum computer, dividing 
the time into small intervals and approximating the time evolution in each slice. The usual procedure utilizes the 
finite n approximation of the Trotter-Suzuki formula and the resulting unitary transformations are applied using a 
circuit-based quantum computer. Second, shows that any quantum computation can be simulated by AQC with 
Tinit polynomial in L (the number of gates in the circuit U). For this, they utilize a Feynman-type Hamiltonian 
similar to ([3]), with a particular encoding of the clock register. The key element is a lower bound on the eigenvalue 
gap of the Hamiltonian (|39p for all t - an inverse polynomial in L. 

However, in j30| . Lloyd has shown that a fast passage from Hinu to the final Hamiltonian works just as well. Let us 
consider the extreme case - an almost instantaneous change in the Hamiltonian. The state of the system then stays 
close to the initial state. The system then continues to evolve according to the final Hamiltonian. As we have seen 
in Section IIII Al this Hamiltonian evolves the initial state within a particular computational subspace (j33p spanned 
by plane- wave eigenvectors of the Hamiltonian of the quantum walk on the line {|^t) ,t = 0, . . . ,L}. Within this 
subspace, it rapidly 'mixes', in the sense that at a random time between and 0(L log L), the system is likely to have 
an active site in the clock register in the success region (see Figure [1]). The sectors of the Hilbert space containing 
IVo)u, ^ \0)c (computing on the proper initial state) and |'Po^)«; "Ki |0)^ (computing on badly initialized work qubits) 
do not mix, because the Hamiltonian does not couple them. Running the HQC then can give us the resulting state 
U \ ipo) with high probability, even though we turned on _ff35 quickly. The proof of AQC's universality thus works just 
as well without the adiabatic theorem (requiring some extra time evolution after reaching the final Hamiltonian). The 
essential element required for this result is that the system evolves only within a preferred computational subspace. 
This in turn is the result of our railroad-switch clock register construction - which is frustration free, and doesn't 



B. HQC and Adiabatic Quantum Computation 



n 




(38) 




(39) 
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induce any 'illegal' transitions. Therefore, our result means that previous universality results for 3-local Hamiltonians 
which required large penalty terms can be greatly simplified, and are valid even without assuming adiabaticity. We 
can now ask where the real power of Adiabatic Quantum Computing lies - e.g. whether the model can stay universal 
if we only use restricted terms in the Hamiltonian as in (3^ . 

C. Protecting the Computation 

The Hamiltonian quantum computation model we just presented is based on two premises. First, we initialize the 
system in a computational basis product state. Second, the Hamiltonian of the system is such that the Schrodinger 
time evolution of this state runs within a particular subspace. Transitions bringing the state of the system out of 
this computational subspace are not present in the Hamiltonian. This allowed us to understand the dynamics of this 
model and bound its required running time. We now ask what are the implications of noise on the dynamics of the 
system and how could we possibly protect the system from moving away from the computational subspace (|33p . 

Let us first look at the structure of the Hilbert space in the light of our Hamiltonian. According to the number of 
active sites in the clock register, the Hilbert space is a direct sum TCq ® Hi © 7^2 ® '^3 © • • • . The orthogonal subspaces 
Hj are not coupled by our Hamiltonian. We divide them into three cathegories. 

First, we would like to avoid the subspaces 7ij>i with more than one active clock site. We can deal with them by 
a 2-local, 0{L) norm, clock-checking Hamiltonian 

i/docfc = E|ll)<ll|c,.c.- (40) 

Its terms add an energy penalty for two neighboring (also within individual railroad switches) active sites in the 
clock. Although states such as l^p)^ ^ |010 • • • 010)^ are not directly violating (PD|) . eigenstates of H containing such 
clock registers have nonzero energy. An eigenstate with two trains in the clock register misses required terms coming 
from the propagation Hamiltonian. If they were all present, terms like \ip')j^ (g) |0 • • • 0110 • • • 0)^ with an illegal clock 
configuration detected by (UOl) would appear. This line of reasoning goes along the lines of the clairvoyance lemma of 
Aharonov et. al. in [13]. It can be used to show that after adding Hdock, the energy of the states in the sectors Hjyi 
is at least 0{L^^). 

Second, we have to deal with Ha - the computationally "dead" subspace with no active site in the clock register. 
This is a problem of the pulse clock (see Section l!!]) . Noise which could kill the active site in the clock could thus 
also kill the computation. One way to prevent this is to energetically favor a single active site by adding a carefully 
balanced Hamiltonian of the form 

Hactive = -ai J2 11) (l|c, + «2 5] |ll) (HI,^,,, , (41) 

where the sum in the second term is taken over neighboring qubits in the clock register. Together with the clairvoyance 
lemma in action, this Hamiltonian favors a single active site over none (or over too many of them). However, adding 
Hactive mcans that our total Hamiltonian can no longer be rewritten as a sum of 3-local projector terms with a 
common zero-energy ground state (and thus be frustration-free) [33| . 

Finally, we have the subspace Hi with a single train, which can be further decomposed into subspaces corresponding 
to different initial states \^po)^ of the n work qubits: 

2"-l 

Til = K^"^. (42) 

Out of these, H^^^ = Hiegai is our computational subspace, corresponding to all states that can be reached by 
Schrodinger time evolution from the initial state (|571) . The eigenvectors of H in the subspace Hi are the momentum 
states \p, ipo), labeled by momentum p and the initial state of the work qubits \^o)w- ^'-'^ badly initialized states 
with ipo 0, Lloyd [30| showed that {p,(pa ^ 0\Hdatainit \p,^a 7^ 0) > 0{L~^). Therefore, the energy gap between 
the computational subspace H^^'' and the subspaces '}{^°^^'> is again on the order of 0{L~^). 

IV. THE SECOND MODEL: A UNIVERSAL 2-LOCAL QUBIT-QUTRIT HAMILTONIAN 

Could we modify the Quantum-3-SAT Hamiltonian of Section IIIII to make it only 2-local? This is probably too 
much to ask when restricting ourselves to qubits. Nevertheless, by using qutrits in place of some of the qubits, we can 
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FIG. 4: The 2-local railroad switch gadget made from qubits (squares) and qutrits (triangles). Simple transitions of the train 
are denoted by black lines. Transitions between two active states of a qutrit are denoted by half-circles. 



obtain a 2-local construction. First, we modify the model from Section IIIH and show that the result is a quantum 
walk on a necklace - a cycle with regularly spaced extra nodes connected to it. Using the results of [s^, we bound 
the required running time of this model. 

We modify the gadget from Figure [3] by using three qutrits and two qubits on each track, as in Figure [H A qubit 
can be active |1) or inactive |0). On the other hand, each qutrit will have two active states |L), |R) (left and right) 
and one inactive state |0). We label the possible positions of the active site in the clock register (the train) on the 
upper track by u\,u\,U2,u)^,u^,U4,u^, and denote the corresponding states of the clock register \u\)c, \u\)c, etc. 
Similarly, we label the possible train positions on the bottom track d\, . . . and the corresponding clock register 
states \l\)c, \li)c, etc. 

There are two types of 2-local transition rules for the upper track. First, we have 6 straightforward rules corre- 
sponding to the thick black lines on the upper track in Figure [4| 

\u\){ct-l\ + \ct-l){u\l \U2){U\\ + \U\){U21 

\u]i){u2\ + \u2){ul\, \u4){ul\ + \ul){u4\, (43) 

|W5>(W4| + |U4)K|, \ct){ul\ + \ul){ct\. 

Second, we have the transitions denoted by red half-circles in Figured They are also 2-local, involving a work qubit 
and a clock qutrit: 

\l){lU,®{\ul){u\\ + \u\){ul\), (44) 

<.®(K)(«3l + l"3>("!l), (45) 

\l){l\,,®{\ul){ul\ + \ul){ul\). (46) 

The rules (|44|) and ((46|) are conditioned on the state of the control qubit qi, and (|45|) flips the target work qubit q2 
when transitioning between jug) and \u^) in the clock register. 

Similar transition rules apply to the bottom track. The straightforward terms are 

\i\){ct^,\ + \c,^,){i\i \i2){n\ + \n){i2i 

\l]i){l2\ + \l2){l]i\, + (47) 

\l^,}{h\ + \h){lh\, \ct){ll\ + \ll){ct\. 

and the analogues of (|^^ and (HS)) for the lower track are conditioned on the control qubit qi being |0), i.e. 

\0){OU,®{ri){l\\ + \l\){ll\), (48) 

|o>(o|,,®(|;^)(;^|-i-|z|)(z||). (49) 



Finally, the analogue of (PSj) is simply 



\ul){ul\ + \ul){ull (50) 



leaving the work qubit q2 untouched. Summing all the terms in (I43p - (|50p gives us the 2-local qubit-qutrit railroad 



switch Hamiltonian H. 



(t) 

23 



Similarly to what we did in Section IIII Al let us examine where the transitions in this Hamiltonian take the state 

l^t^i) - (a |0)„, +b |1)„, ) ® |ci_i), (51) 
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FIG. 5: The geometry of connections in ffjs restricted to the subspace 'H„23 (within one gadget) 




FIG. 6: The necklace graph depicting the connections in H23 in a particular basis. The shaded region corresponds to states 
within one railroad-switch gadget. 



where a |<^o...) + b I'fii...) is a general state of the work register and the active site of clock register is at the qubit Ct-i- 
Within the gadget, the transitions in H23 move this state around a 12-dimensional subspace 7ig*32 with the geometry 
of allowed transitions depicted in Figure [5l The basis of this subspace contains first the endpoint states 



|C\> - {a\^o...)+b\^i...))®\ct-i] 



where \ip'i ) = a^Jipi,,} has the target qubit (72 flipped. Second, we have the sequence 



(52) 
(53) 

(54) 
(55) 
(56) 
(57) 
(58) 
(59) 
(60) 
(61) 

Finally, the last two states in the basis are the two 'blind-alley' states, in which the active site tries to use the track 
which it is not supposed to take 





= a 


Vo...) 






® hi),, 




— a 


Vo...) 








m) 


— a 


^0...) 


1^2),- 


^b\ipi...) 






= a 


^0...) 


® 1^3),- 


\^b\(pi,„) 






= a 


^0...) 




^M^'i...) 






— a 


Vo...) 


® 1^4), ' 


^b\ip[,„) 






— a 


^0...) 




^b\ip[,„) 






= a 


^0...) 




\-b\ip[,„) 





1*), 

-1 /s 



a\(po...)<»\u\)^ + b\<fi,„)<»\l\)^, 
a\^o...)®\ul)^ + b\ip[J®\ll)^. 



(62) 
(63) 



Omitting the connection to the region outside the gadget, this 12-dimensional subspace Hg*23 is invariant under the 



transitions in i?23^- is convenient to analyze the dynamics of H2^' using the basis just given. In this basis, the 



(t) 



Hamiltonian within one gadget has a very simple form - the adjacency matrix of a line , \ip. 



''t-i/, ■ 



with two additional nodes {ipt^i) and I'f/'j^*/) as in Figure [5] Moreover, we can give our clock register a very regular 
form by arranging 8 clock qubits between two successive gadgets. We can then choose to look at the dynamics of 
our system in a basis consisting of states of the form ^ (for the non-gadget parts) and of the bases just given for 
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individual gadgets. Similarly to what we found in Section [ill Al in this basis the Hamiltonian of our model has a very 
simple form. It is an adjacency matrix of a ring with extra lines sticking out of it on every 9-th point as in Figure [S] 
This allows us to forget about the content of the work register when solving the dynamics of the time evolution with 
our Hamiltonian. 

Finally, we need to find the required running time for the qubit-qutrit model. The dynamics of the time evolution 
generated by the Hamiltonian is a quantum walk on a necklace from Figure [6l In an upcoming paper by Nagaj and 
Reitzner I34|, we investigate the mixing properties of this type of quantum walks. Again, the mixing is fast. When 
we start in the initial state and let the system evolve with iJ23 for a random time < t < T with T — 0{L^) 
(up to logarithmic factors) and measure the output work qubit, our model is universal for BQP. 
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APPENDIX A: CONTINUOUS-TIME QUANTUM WALKS IN ID 

First, in Section [A II we analyze the continuous-time quantum walk on a line and prove a lemmas about the mixing 
of this walk. Second, in Section [A 21 we analyze the quantum walk on a cycle and prove another mixing lemma used 
in the proof of universality of the railroad switch Hamiltonian Computer in Section IIIII 



1. Quantum Walk on a Line 



Consider a continuous time quantum walk on a line of length L, where the Hamiltonian is the negative of the 
adjacency matrix for the line 

i?i = -E + + + (Al) 



The eigenvalues of this Hamiltonian are 

A, — —2 cos I 

for j — 1 . . . while the corresponding eigenvectors |(/)^-'-') — \k) have components 



^J=-2cos(-^), (A2) 



Consider the time evolution of a particular basis state |c). The probability of finding the system in the basis state 
|m) at some time r can be found by expanding |c) and \m) in the basis of the eigenvectors (jASP : 

L 

Pr{m\c) = I (ml e-^^- |c)f - ^ e-''(^^-^'=)-0(^)0(^>0W>W . (A4) 

Because the time evolution (according to the Schrodinger equation) is unitary, this probability pr{m\c) does not 
converge. On the other hand, let us define the time average of pr{m\c) for time < r < T20 as 

Pr2o(jn\c) = — / Pr{m\c)dT. (A5) 
T20 Jo 
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As we will show below in Lemma [1] this average probability distribution does converge to a limiting distribution 
7r(m|c), defined as the T20 00 limit of the average probability distribution (|A5D . All the eigenvalues (IA2p are 
different, so we can express the limiting distribution as 



L 

^(m|c)= hm P..o(m|c)=^|0Hf |#r, (A6) 



which in this case is 



-He) = . (A7) 



According to the following lemma, the average probability distribution (|A5[) converges to the limiting distribution 
7r(m|c). 

Lemma 1. Consider a continuous time quantum walk on a line of length L, where the Hamiltonian is the negative of 
the adjacency matrix for the line. Let the system evolve for time r < chosen uniformly at random, starting in a 
position basis state |c). The average probability distribution Pt2o('|c) converges to the limiting probability distribution 
7r(-|c) as 

J2 \PrMc) - nim\c)\ < O (—) . (A8) 



ni—l 



T20 



Proof. First, recall Lemma 4.3 of pBa ] for the total variation distance of the probability distribution p^-^^ from the 
limiting distribution, saying 

^ \prMc) - n{m\c)\ < — E TT^ri- (^9) 

Using (jA2p and (jA3p . we can bound the expression on the right of (jA9p . When j is close to k, i.e. |j — A:| < Ci, we 
can write 



a(j)|2 



|Aj - Afcl 

On the other hand, for \ j — k\ > Ci we can bound 



< 2. (AlO) 



< (AH) 



\\j~Xk\ L + V 
with Ci and C2 constants independent of L. Inserting into (jA9[) . we have 

E \PrMc) - Am\e)\ < ^ + ^ = O f A) . (A12) 

m— 1 ^ 

□ 

2. Quantum Walk on a Circle 

If the geometry of the system is a closed loop of length L instead of a line, the Hamiltonian (|A1[) gets an additional 
wrap-around term. 

L-l 

H,o,, = - {\L) (1| + |1> {L\) - E (b) U + 1| + \J + 1> 01) ■ (A13) 
The eigenvalues of this Hamiltonian are 

Aj = -2cospj, (A14) 
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corresponding to plain waves with momenta 



for j = ... L — I. The corresponding eigenvectors 1'/)^-''') = J2k=i 'rk 1^) have components 



= (A15) 



V ^ 



^fc - A/T^'^^ ' J = (A17) 



These can be combined to make real eigenvectors. For our analysis, it will be enough to consider a line with even 
length, and only the cosine plain waves: 



VI' 

ij) [2 , . _ L 



(A18) 

k - Y;^cos(pjfc), j = l,...,-^. (A19) 

The limiting distribution on the cycle when starting from site c is 

1 2 

Z " Z 

for all points m except for m = c (return back) and m = L + 1 — c (the point across the cycle), where we have 

7r(m|c) = I - (A21) 

As in Section [A H we will again utilize Lemma 4.3 of 35] to prove the convergence of the time-averaged probability 
distribution to this limiting distribution. On the right side of (jA9p . we now have 



7r(m|c) = - - — . (A20) 



|4^T<|- (A22) 

The sum over the non-equal eigenvalues 

^° ^ E TT^i (A23) 

is now more complicated, because of the degeneracy of the spectrum. We plot |Aj — Afc|~^ in Figure [T] Because of the 
symmetrical way Xj (|A14[) arise, assuming L is divisible by 4, we can express (|A23|) as 

L/4 L/2-J 

j=0 k=j+l I J *= I 

L/2 L/4L/2-j 

= ^^y^^ ^+16E E R (A26) 

' «- ' ^ « ' 

The term Aq comes from j = 0. We can bound the sum by an integral, taking x = obtaining 



A, 



it^ 1 L dx 







- 1 - cos (27r^) ~ 2" 1 - cos(27ra;) ' ^ ^ 







L i 



FIG. 7: A contour plot of \Xj - Xk\ ^ with \j from (|XT4)) . The sum in (|A25P is over the marked region. 



For X £ [O, ^] , we can bound 
resulting in 



1 — cos 27rx > 2ttx 



We bound the other term, Ai, in (|A26p by an integral as well: 

L/4 L/2-j 



^Uk^+i cos(27rA)-cos(27ri) 
" dx 



< — / rfy 



cos (27ra;) — cos (27r?/) 



Again, we can lower bound the difference in eigenvalues for y G [O, ^] and x £ [v^ h ^ 

cos {2nx) - cos {2ny) > 2it{x^ - y"^), 



allowing us to write 



L2 
47r 

47r 



1 _i 

4 Z, 



L2 /-I-T 
477 



2y 



9 9 
+y X - 



2y 



[log(l-4y) + log(l + 2yi)]. 



As y in (IA35P is at most 
Finally, this results in 



4 L 



we can bound R by 



|i?| <logL. 



L^logX /-"T L^logL , 

< / — = log 

27r ii 2?/ 27r ^ 



1 _ j_ 

4 L \ ^ ^^r2 1„„2 



L 

<logL 



< 0(i^ log' L). 



(A28) 
(A29) 

(A30) 
(A31) 

(A32) 

(A33) 
(A34) 
(A35) 

(A36) 
(A37) 
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Putting (K29\i and (|M7| into (|M6)) . we obtain 



Lemma 4.3 of [Hj (see (|A9|) ) then reads 



S-o < 0{L^\og^L). 



2 1 



m|c) - 7r(m|c)| < — -0(L^ log^ L) = O 

To L 



i log^ L 



(A38) 



(A39) 



Thus, for Tcirc = eO(Lfog^ L), the time-averaged distribution becomes e-close to the hmiting distribution. Using the 
bound on the total variation distance we just proved, it is straightforward to obtain the foUowing lemma which we 
use m Section [ml 

Lemma 2. Consider a continuous time quantum walk on a cycle of length L (divisible by 4), where the Hamiltonian 
is the negative of the adjacency matrix for the cycle. Let the system evolve for a time t < Tq chosen uniformly at 
random, starting in a position basis state |0). The probability to measure a position state \t) farther than L/6 from 

the starting point (the farther two thirds of the cycle) is then bounded from below as po > ^ — — O ^ ^ ^"g ^ ^ _ 

Proof. We proceed as in the proof of Lemma[T]in Section IXTl Let us call the farther two thirds of the cycle (see Figure 
[T]in Section HITI the success region (SR). When we choose the time t < Tq uniformly at random, the probability to 
measure a state \t) with t e SR is 



(A40) 



meSR 



Using the bound on the total variation distance (|A39[) we just proved and the formulae for the limiting distribution 
(|MOl) .(|MI |) . we have 



O 



L loK^ L 



?7l — 1 

> Y IPto(to|c) -7r(TO|c) 



m£SR 



> 



meSR 



meSR 



3L 



L 



(A41) 
(A42) 

(A43) 
(A44) 



Therefore, the probability of finding the chain in state IV'teSi?) at a random time r < Tq is bounded from below by 



2 1 

^°^3-3Z-^ 



L log" L 



(A45) 



□ 



It is thus enough to wait a random time not larger than 0(Llog^ L) to find the state of the system in the success 
region with probability close to two thirds. 
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